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ABSTRACT 

The  purpose  of  this  report  is  to  present  graphical  solutions  to  the 
equations  for  isentropic  flow  of  ideal  gases  at  subsonic,  sonic,  and 
supersonic  velocities.  Knowing  only  the  weight  flow,  molecular  weight, 
pressure,  temperature,  k  values,  and  the  exit  pressures  (vacuum),  the 
reader  will  be  able  to  design  a  flow  system  by  specifying  pipe  sizes, 
venturi  throats,  and  exit  diameters  without  the  use  of  computers. 
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1.  Introduction 


In  laser  work  optical  energy  is  extracted  from  thermally 
excited  species  of  fast  flowing  hot  gases.  For  this  extraction  a  number 
of  ambient  temperature  gases  must  be  moved  in  ducts,  the  mass  flow  must 
be  controlled  by  venturis,  and  they  must  be  combined  to  burn  at  low 
pressure  and  high  velocity  in  the  cavity. 

Design  calculations  must  be  made  for  every  line  and  every  venturi. 
Despite  the  relative  simplicity  of  the  calculations,  they  are  time 
consuming.  A  nomographic  procedure  enables  engineers  to  make  fast  cal¬ 
culations  without  use  of  a  computer.  The  graphs  in  the  report  are 
designed  for  conmonly  used  variables  such  as  temperatures  in  degrees 
Fahrenheit. 


2.  Discussion 

It  is  assumed  that  the  flow  in  the  ducts  is  adabatic.  The 
basic  equations  to  describe  the  flow  are  taken  from  "Compressible  Fluid 
Flow,"  Volume  I,  by  A.  H.  Shapiro,  and  are  as  follows. 


(1) 

4.16) 


=  (l  +  M2^  1  (p  83,  Eq.  4.14b)  (2) 

co  * k  *  E  To(sr^)  <p  79- E"-  4-5b)  (3) 
C*2  *  «  kTT  R  To(moTut  )  <»  80’  E<>'  4-5,::)  <4> 

V  =  M*  C*  (p  81)  (5) 


1 


The  graphical  solutions  of  these  equations  (presented  in  the  report) 
are  accurate  but  inaccuracies  using  a  straight  edge,  marking  the  refer¬ 
ence  lines  and  going  to  the  next  line,  will  degrade  them  so  that  the 
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solutions  will  be  only  within  ±10  percent.  The  advantage  of  this  method 
of  calculation  is  that  it  can  be  done  without  the  use  of  mathematics 
and  that  it  is  almost  impossible  to  make  a  mistake  by  a  factor  of  10  or 
100  which  is  easy  to  do  on  a  slide  rule.  A  further  advantage  is  that 
trends  can  very  easily  be  observed  by  varying  certain  key  inputs  within 
limits . 

The  equations  are  correct  for  ideal  gases  only,  but  they  will  give 
acceptable  results  for  gases  at  low  pressure  and  at  temperatures  high 
above  the  critical  temperature.  They  are  also  good  for  mixtures  of 
gases  (such  as  air)  if  the  average  specific  heat  ratios  k,  and  the 
average  molecular  weights  are  used. 

If  burning  takes  place,  there  is  a  different  set  of  gases  before  and 
after  the  reaction.  The  temperature  will  change  due  to  the  heat  of 
reaction,  the  number  of  moles  will  not  be  the  same,  and  the  value  of  the 
specific  ht  't  ratio  will  also  change. 

For  example  the  burning  of  carbon  disulfide  with  oxygen  to  carbon 
monoxide  and  sulphur  dioxide  uses  7  moles  of  reactants  to  produce  6  moles 
of  products  as  shown: 

2CS2  +  5C2  =  2C0  +  4S02 

The  k  values  for  the  chemicals  are  as  follows: 

CS£  -  1.21 

CO  -  1.40 
02  -  1.40 

S02  -  1.29 

The  molecular  weignt  of  the  products  is  7/6  times  that  of  the  reactants. 

Furthermore  the  heat  ratios  change  with  temperature  and  pressure 
and  most  experimental  data  are  only  for  ambient  temperatures.  The 
International  Critical  Tables,  Volume  5,  pp.  80-82,  show  for  example, 
that  for  oxygen  at  1  atmosphere  the  k  value  increases  from  1.40  to  1.45 
as  the  temperature  drops  from  15°C  to  -181 °C.  Using  Graph  3  to  calculate 
the  pressure  at  Mach  1,  the  error  introduced  would  be  approximately 
1.5  percent.  However,  if  experimental  data  for  the  reaction  products  at 
elevated  temperatures  are  available  the  equations  may  be  solved  for  the 
flow  problems. 
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3.  Description  o^  the  Nomographs 

The  variables  in  Graph  1  are  k  and  M. 
factor  (1  +  (k  -  1)  M  •  M/2)  common  to  Graphs  2, 


The  solution  is  the 
3,  and  6. 
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Graph  2  is  the  solution  of  Equation  (1).  Straightforward  use  solves 
for  the  flow  area  but  any  other  variable  may  be  found  as  shown  in  the 
following  example.  For  flow  velocities  slower  than  Mach  0.4,  it  is  not 
necessary  to  solve  for  the  factor  in  Graph  1.  The  calibration  on  the 
lower  part  of  line  2  is  sufficient  for  the  solution  but  accuracy  is 
impaired.  Graph  3  solves  Equation  (2);  Graph  4  solves  Equations  (3), 
(4),  and  (5);  and  Graph  6  solves  Equation  (6). 
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Graph  2  can  be  used  straightforward  if  the  unknown  is  the  area.  ■* 
If  the  unknown  is  any  other  variable,  the  procedure  is  the  same  until  the  | 
unknown  is  reached  then  the  solution  is  exactly  backwards  starting  with  J 
the  area.  For  low  velocities  (less  than  Mach  0.4),  line  2  is  calibrated  | 
in  Mach  numbers  and  this  should  be  the  starting  pcint  rather  than  line  A.  | 


Graph  3  solves  for  the  flowing  pressure  and  is  more  accurate  at 
velocities  greater  than  Mach  1.  Note  that  at  the  intersection  of  line  B 
and  line  1  the  total  pressure  and  the  flowing  pressure  are  equal.  If  the 
pressure  range  is  not  on  this  graph  (like  3100  mm  Hg)  use  31  and  multiply 
the  answer  by  a  factor  of  100  as  well  (read  16.4  x  100  =  1640  mm  Hg  for 
air  at  Mach  1.0). 

Graph  4  solves  for  the  gas  velocity  at  Mach  1  (in  ft/sec)  and  then 
by  straight  multiplication  for  any  other  velocity.  The  conversion  chart 
shows  the  distance  traveled  in  mm/^sec. 

Graph  5  shows  a  comparison  of  inlet  and  exit  diameters  to  the 
throat  diameters  for  k  values  ranging  from  1.1  to  1.66.  It  is  a  solution 
of  Equation  1. 


Graph  6  gives  the  temperature  at  the  flow  velocity  in  conjunction 
with  Graph  1. 

The  development  of  the  graphs  using  the  formulas  indicated  was 
straightforward.  The  graphs  were  prepared  using  the  technique  outlined 
in  "Easy  Way  to  Make  Nomographs,"  Chemical  Engineering,  November  1952. 

Instead  of  calibrating  the  scales  with  their  actual  values,  input 
values  such  as  k  were  used.  For  example  in  Graph  1  the  specific  heat 
ratio  scale  is  calibrated  from  k  **  1.1  to  k  =  1.65.  The  actual  value 
at  the  point  1.1  is  0.05  which  is  (1.1  -  l)/2  and  at  k  =  1.6  it  is  0.3. 
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The  slanting  line  in  Graphs  2  and  3  occurs  because  of  raising  the 
result  from  Graph  1  to  a  power.  Again  the  k  scale  (slanting  line)  is 
calibrated  in  terms  of  k  rather  than  (k  +  l)/2(k  -  1)  as  in  Graph  2, 
or  k/(k  -  1)  as  in  Graph  3. 

The  temperature  scales  in  Graphs  2  and  4  are  calibrated  in  °F.  The 
actual  value  at  0°F  is  460  (°R)  and  at  2000°F  it  is  2460  (°R).  The 
length  of  1  log  scale  is  from  100°R  (-360°F)  to  1000°R  (540°F). 


4.  Method  of  Solution 

Sample  solutions  using  each  and  every  graph  will  show  how  to 
solve  problems. 

Air  was  used  in  the  problem  because  the  reference  book  furnishes 
short-cut  formulas  for  solving  the  equations.  A  numerical  solution 
follows  the  graphical  one  so  that  the  reader  may  compare  results.  It 
should  be  noted,  however,  that  the  graphs  will  solve  problems  for  any 
gas  while  the  short-cut  formulas  are  valid  only  for  air. 

EXAMPLE: 

The  following  data  are  given: 


Working  fluid 

Air 

Air 

Molecular  weight 

f 

28.9  g/mol 

28.9  lb/ lb  mol 

Weight  flow  (W) 

45.4  g/sec 

0.1  lb /sec 

Stagnation  temperature 

1140°F 

1600°R 

Stagnation  pressure 

3100  mm  Hg 

60  psia 

Specific  heat  ratio  (k) 

1.4 

1.4 

Universal  gas  constant  (R) 

1544 

Gravitational  constant  (g) 

32.2  ft/sec2 

The  following  are  required: 


a)  The  throat  diameter  of  a  critical  flow  venturi 

b)  The  pressure  at  the  throat 

c)  The  temperature  at  the  throat 

d)  The  velocity  of  sound  at  the  stagnation  conditions 


e)  The  actual  velocity  of  the  gas  flowing  through  the  throat 

f)  The  velocity  of  the  gas  in  a  1-inch  diameter  (ID)  line. 


a.  Solution 


1)  At  a  sonic  throat  the  velocity  is  always  Mach  1.0. 

On  Graph  1  connect  k  =  1.4  with  M  =  1.0  and  read  1.2.  On  Graph  2 
connect  1.2  on  line  A  with  1.4  on  line  B  and  mark  line  1.  Connect  the 
point  on  line  1  with  the  Mach  number  1.0  on  line  A  (note  that  line  A  is 
used  twice,  and  for  different  values  each  time)  and  mark  line  2.  Connect 
with  the  pressure  of  3100  mm  Hg  on  line  C  and  mark  line  3.  Connect  with 
the  weight  flow  of  45.4  g/sec  on  line  D  and  mark  line  4.  Connect  with 
the  specific  heat  ratio  k  =  1.4  on  line  E  and  mark  line  5.  Connect  with 
the  temperature  of  1140°F  on  line  F  and  mark  line  6.  Connect  with  the 
molecular  weight  on  line  G  and  read  the  flow  area  as  0.13  square  inch. 

The  horizontal  line  below  gives  the  throat  diameter  as  0.41  inch. 

2)  For  the  pressure  at  the  throat,  connect  on  Graph  1 
k  =  1.4  with  M  =  1.0  and  read  1.2.  On  Graph  3,  connect  1.2  on  line  A 
with  1.4  on  line  B  and  mark  line  1.  Connect  with  line  marked  "Total 
Pressure"  (3100  mm  Hg)  and  read  1640  mm  Hg  on  the  line  marked  "Pressure." 

3)  For  the  temperature  at  the  throat,  use  again  the 
factor  1.2  from  Graph  1.  On  Graph  6,  connect  this  number  with  1140°F 
on  the  line  marked  "Stagnation  Temperature"  and  read  890 °F  on  the  line 
marked  "Flowing  Temperature." 

4)  The  velocity  of  sound  is  calculated  by  Equation  (3) 
and  is  the  velocity  of  the  propagation  of  sound  in  a  medium  at  rest. 

On  Graph  4  the  left-hand  scale  of  line  A  (k  to  calculate  Cq)  is  used, 

connect  1.4  with  28.9  on  line  B  and  mark  line  C.  Connect  with  1140°F 
on  line  D  and  read  1950  ft/sec  on  line  E. 

5)  The  actual  velocity  of  the  gas  flowing  through  the 
throat  is  given  by  the  same  graph  (Graph  4)  but  using  the  right-hand 
scale  on  line  A  (k  to  calculate  C*).  Connect  1.4  on  line  A  (right-hand 
scale)  with  28.9  on  line  B  and  mark  line  C.  Connect  with  1140°F  on  line 
D  and  read  1750  ft/sec  velocity  on  line  E. 

6)  The  velocity  of  the  gas  in  a  1-inch  ID  line  is  sub¬ 
sonic  because  the  throat  diameter  (for  sonic  flow)  is  0.41  inch.  Because 
the  flow  diameter  and  thereby  the  area  are  given  and  the  velocity  of  the 
gas  is  unknown.  Graph  2  must  be  solved  backwards.  On  the  horizontal 

2 

line  1-inch  diameter  corresponds  to  0.785-inch  area.  Connect  0.785  on 
line  marked  "Flow  Area"  with  28.9  on  line  G  and  mark  line  6.  Connect 
with  1140°F  on  line  F  and  mark  line  5.  Continue  until  line  3  is  marked. 
Connect  with  3100  inn  Hg  pressure  on  line  C  and  read  Mach  0.09  on  line  2. 
The  value  of  C*  (1750  ft/sec)  should  be  multiplied  by  0.09  to  give  a 
velocity  of  157  ft/sec  on  the  1-inch  line. 
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b.  Numerical  Solution 

1)  For  the  throat  diameter,  Flieguers  formula 
(p  85  -  4.18)  is  used  with  the  English  units  (such  as  W  =  0.1  lb/sec) 
called  for  in  the  text. 


V*0 

PQ  X  0.532 


.1  lb/sec S/1600°R 
60  psia  (0.532) 


=  0.125  in.‘ 


q2  =  0.125  x  4 

JT 


0.1585 


D  =  0.4  in. 

2)  For  the  pressure  of  the  throat,  the  formula 
P*/Pq  =  0.5283  can  be  used  in  nm  Hg  absolute  because  P*/P_  is  a  ratio 

P*  =  3100  x  0.5283  =  1635  mm  Hg. 

3)  The  temperature  at  the  throat  is  again  in  ratio  form 
T*/Tq  =  0.8333,  but  absolute  units  such  as  Rankin  or  Kelvin  must  be  used. 

T*  =  T(0.8333)  =  (1600) (0.8333)  =  1335°R 

T*  =  1335  -  460  =  875 °F 

4)  The  velocity  of  sound  is  given  by 

C?  =  k  g  R  Tn  7  .1  -r 

0  °  0  (mol  wt) 

1 

(  1.4  x  32.2  x  1544  x  160o\2  , 

CQ  - 28T§ - )  =  1963  ft/s6C  * 
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.Vtul  iil.l  >•'*' : 


A  =  0.0726 


(l  +  0.2  M2) 

M 


=  0.785 


Since  the  graphical  solution  gave  M  =  0.09,  this  will  be  the  assumed 
starting  point. 


M 

0.2  M2 

(l  +  0.2  M2) 

(l  +0.2  M2) 

M 

3 

0.0726  —  *  °*2  M^- 

| 

0.09 

0.00162 

1.00486 

11.1 

0.805 

f 

0.095 

0.00182 

1.00547 

10.55 

0.77 

J 

i 

0.100 

0.002 

1.006 

10.06 

0.73 

1 

|  1 

,j  1 


Since  the  actual  flow  area  is  0.785  in.  ,  the  Mach  number  must  be 
between  0.09  and  0.095. 

Assuming  M  -  0.093, 


A  =  0.0726  x 


[l  +  (0.2)(0.093) 
(0.093) 


=  0.785 


The  velocity  then  is  M*  C*  =  0.093  x  1792  =  167  ft/sec. 


rjV  I  ^ J  -57  y  ^!'A?  j  jrigVjiX;^  W^H<r^^A- W1-  ?  f  V,'Tii^^^>7  r>^  * 


I 

> 


f 


: 

t 


v 

p 

I 


\  3.0  4.0 

■ii, ni  l..  1. 1  i  ii  1. 1 1  li  i  1 1 1 

H  |  '  1  '  1  I 

t-  7.0  10.0 


10,000  °R 


GRAPH  3 


VALUE 


I 


f  Si,’»6AJI 
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CALCULATION 
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THIS  GRAPH  SOLVES  THE  EQUATION  AND  f 
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v0  =  c0  x  M  WHERE  C§  =  k  g  R  JqI  MOL  WEIGHT  '( 
V*  =  C*  x  M  WHERE  C*2  =  k  g  R  TV  MOL  WEIGHT  ■] 


CONNECT  LINE  A  WITH  B,  MARK  C  (REF  LINE  1).  GO  T t 
D,  READ  AND  MARK  E.  THEN  GO  TO  F  AND  READ  VEl| 


IS  GRAPH  SOLVES  THE  EQUATION  AND 


i  ■  CQ  x  M  WHERE  Cjj  »  k  g  R  Tq/  MOL  WEIGHT 
V*  =  C*  x  M  WHERE  C*2  =  k  g  R  T*/  MOL  WEIGHT  =  R  Tq/MOLES 

lx  >  I 


NNECT  LINE  a  WITH  B,  MARK  C  (REF  LINE  1).  GO  TO 
READ  AND  MARK  E.  THEN  GO  TO  F  AND  READ  VELOCITY. 
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